arXiv:hep-ph/0211292vl 19 Nov 2002 


Bled Workshops 
IN Physics 
VOL. 3, No. 3 


Proceedins of the Mini-Workshop 
Qiiarks and hadrons (p. 67) 
Bled, Slovenia, July 7-14, 2002 


A 


Axial currents in electro-weak pion production 
at threshold and in the A-region * 

S. Sirca“’^, L. Amoreira‘^’‘^, M. Fiolhais^ ®, and B. Golli^’^ 

'^Faculty of Mathematics and Physics, University of Ljubljana, 1000 Ljubljana, Slovenia 
'’Jozef Stefan Institute, 1000 Ljubljana, Slovenia 

■^Department of Physics, University of Beira Interior, 6201-001 Covilha, Portugal 
■^Centre for Computational Physics, University of Coimbra, 3004-516 Coimbra, Portugal 
‘■Department of Physics, University of Coimbra, 3004-516 Coimbra, Portugal 
^Faculty of Education, University of Ljubljana, 1000 Ljubljana, Slovenia 


Abstract. We discuss electro-magnetic and weak production of pions on nucleons and 
show how results of experiments and their interpretation in terms of chiral quark models 
with explicit meson degrees of freedom combine to reveal the ground-state axial form 
factors and axial N-A transition amplitudes. 


1 Introduction 

The study of electro-weak N-A transition amplitudes, together with an under¬ 
standing of the corresponding pion electro-production process at low energies, 
provides information on the structure of the nucleon and its first excited state. 
For example, the electro-magnetic transition amplitudes for the processes 'Y*p —> 
A+ —> p7T° and Y*p —> A+ —> n7T+ are sensitive to the deviation of the nucleon 
shape from spherical symmetry [|l[]. Below the A resonance (and in particular 
close to the pion-production threshold), the reaction Y*p —> n7T+ also yields in¬ 
formation on the nucleon axial and induced pseudo-scalar form-factors. While 
the electro-production of pions at relatively high and low momentum 
transfers has been intensively investigated experimentally in the past years at 
modern electron accelerator facilities, very little data exist on the corresponding 
weak axial processes. 

2 Nucleon axial form-factor 

In a phenomenological approach, the nucleon axial form-factor is one of the quan¬ 
tities needed to extract the weak axial amplitudes in the A region. There are ba¬ 
sically two methods to determine this form-factor. One set of experimental data 
comes from measurements of quasi-elastic (anti)neutrino scattering on protons, 
deuterons, heavier nuclei, and composite targets (see for a comprehensive list 
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of references). In the quasi-elastic picture of (anfi)neutrino-nucleus scatfering, the 
"vN —> pN weak transition amplitude can be expressed in terms of the nucleon 
electro-magnetic form-factors and the axial form factor Ga- The axial form-factor 
is extracted by fitting the -dependence of fhe (anti)neutrino-nucleon cross sec¬ 
tion, 

^ = A(Q2] t B(Q2) (s - u) + C(Q2) (s - u)2 , (1) 

in which Ga(Q^) is contained in the A(Q^), B(Q^), and C(Q^) coefficients and is 
assumed to be the only rmknown quantity. It can be parameterised in terms of an 
'axial mass' Ma as 

GA(Q") = GA(0)/(1+QVMi)^ 

Another body of data comes from charged pion elecfro-production on pro¬ 
tons (see [Q] and references therein) slightly above the pion production thresh¬ 
old. As opposed to neutrino scattering, which is described by the Cabibbo-mixed 
V — A theory, the extraction of fhe axial form facfor from elecfro-production re¬ 
quires a more involved theoretical picture [^H]. The presently available most pre¬ 
cise determination for Ma from pion elecfro-production is 

Ma = (1.077 ±0.039) GeV (2) 

which is AMa = (0.051 ± 0.044) GeV larger than the axial mass Ma = (1.026 ± 
0.021) GeV known from neutrino scattering experiments. The weighted world- 
average estimate from electro-production data is Ma = (1 -069 ± 0.016) GeV, with 
an excess of AMa = (0.043 ± 0.026) GeV with respect to the weak probe. The 
~ 5 % difference in Ma can apparently be attributed to pion-loop corrections to 
the electro-production process [^]. 


3 N-A weak axial amplitudes 


The experiments using neutrino scattering on deuterium or hydrogen in the A 
region have been performed at Argorme, GERN, and Brookhaven WMB- 
(Additional experimental results exist in the quasi-elastic regime, from which Ma 
has been extracted.) For pure A production, the matrix element has the familiar 
form 

M = (pA|vN) = ^^^^j«(A|V“-A“|N), 

where Gp is the Fermi's coupling constant, 0c is the V^d element of the CKM 
matrix, )„ = u^Ya)! — Y 5 )u^ is the matrix element of fhe leptonic current, and 
the matrix element of fhe hadronic currenf has been split into its vector and 
axial parts. Typically either the A++ or the A+ are excited in the process. The 
hadronic part for the latter can be expanded in terms of weak vector and axial 
form-factors [|l^] 


M = ■^ua«(p') 


pV 

'“3 


pv 

'“4 
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where — q“j^^, ua(x(p^) is the Rarita-Schwinger spinor describing 

the A state with four-vector p', and u(p) is the Dirac spinor for the (target) nu¬ 
cleon of mass M with four-vector p. (In the case of the A++ excitation, the expres¬ 
sion on the RHS acquires an additional isospin factor of Vi since (A++ | | p) = 

>/3(A+ I J“ I p) = •\/3(A° 11“ I p).) The function f (W) represents a Breit-Wigner de¬ 
pendence on the invariant mass W of the Ntt system. 

The matrix element is assumed to be invariant under time reversal, hence all 
form-factors Cy’^(Q^) are real. Usually the conserved vector current hypothesis 
(CVC) is also assumed to hold. The CVC cormects the matrix elements of the 
strangeness-conserving hadronic weak vector current to the isovector component 
of the electro-magnetic current: 

(A++|V“|p)=V3(A+|JgM(T=l)|p), 

(A°|V“|p)= (A+|Je“m(T = 1)Ip). 


The information on the weak vector transition form-factors is obtained from 
the analysis of photo- and electro-production multipole amplitudes. For A electro¬ 
excitation, the allowed multipoles are the dominant magnetic dipole Mi+ and 
the electric and coulomb quadrupole amplitudes Ei+ and Si+, which are formd 
to be much smaller than Mi + [^H]. If we assume thaf Mi + dominates the electro¬ 
production amplitude, we have C 5 = Cg = 0 and end up with only one indepen¬ 
dent vector form-factor 


v__M 

4- 


It turns out that electro-production data can be fitted well with a dipole form for 

pV 

'“3' 


C^(q2) ^2.05 



0.54 GeV^ 



An alfernative parameferisation of C 3 which accounts for a small observed devi¬ 
ation from the pure dipole form is 


C^(Q2) 


2.05 


1 +9v^ 


exp 




The main interest therefore lies in the axial part of the hadronic weak current 
which is not well known. 


Extraction of Cf (Q^) from data 

The key assumption in experimental analyses of the axial matrix element is the 
PC AC. It implies that the divergence of the axial current should vanish as —> 

0, which occurs if fhe induced pseudo-scalar term with Cg (the analogue of Gp 
in the nucleon case) is dominated by the pion pole. In consequence, can be 
expressed in terms of the strong ttNA form-factor, 

C^(Q^) /2 G^na 1 

M2 "V 3 2M Q2 m2 ’ 
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while and can be approximately connected through the off-diagonal Gold- 
berger-Treiman relation In a phenomenological analysis, C^(Q^), C 4 (Q^), 
and Cg (Q^) are taken as free parameters and are fitted to the data. The axial 
form-factors are also parameterised in "corrected" dipole forms 


Cf(Q2) =Cf(0) 


uiQ 


2 1 r 


bi + Q' 


1 + 


Mij 


In the simplest approach one takes Ui = bt = 0. Historically, the experimental 
data on weak pion production could be understood well enough in terms of a 
theory developed by Adler For lack of a better choice, Adler's values for 
Ch( 0 ) have conventionally been adopted to fix the fit-parameters at = 0 , i. e. 


c^(0) = 0 , (3) 

C^(0) = -0.3 , (4) 

C^(0) = 1.2. (5) 


In such a situation, one ends up with Ma as the only free fit-parameter. 

Several observables are used to fit the Q^-dependence of the form-factors. 
Most commonly used are the total cross-sections cr(E-j,), and the angular distribu¬ 
tions of the recoiling nucleon 
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where Pmn are the density matrix elements and Y^m are the spherical harmonics. 
Better than from the Pran coefficients, the dependence of the matrix element 
can be determined from the differential cross-section dcr/dQ^. In particular, since 
the dependence on and C 4 is anticipated to be weak at ~ 0 , then 


dcr 

d^ 


(Q" 


0 )oc(Cf( 0 )) 2 . 


The refinements of this crude approach are dictated by several observations. 
If the target is a nucleus (for example, the deuteron which is needed to access 
specific charge charmels), nuclear effects need to be estimated. Another important 
correction arises due to the finite energy width of the A. In addition, the non-zero 
mass of the scattered muon may play a role at low Q^. 

All these effects have been addressed carefully in [p^]. The sensitivity of the 
differential cross-section to different nucleon-nucleon potentials was seen to be 
smaller than 10% even at < 0.1 GeV^. In the range above that value, this 
allows one to interpret inelastic data on the deuteron as if they were data obtained 
on the free nucleon. The effect of non-zero muon mass is even less pronounced: 
it does not exceed 5 % in the region of ~ 0.05 GeV^. The energy dependence 
of the width of the A resonance was observed to have a negligible effect on the 
cross-section. The final value based on the analysis of Argonne data is 

C^(0) = 1.22 ±0.06. 


(6) 
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At present, this is the best estimate for 05 ( 0 ), although a number of phe¬ 
nomenological predictions also exist [p^. We adopt this value for the purpose of 
comparison to our calculations. There is also some scarce, but direct experimental 
evidence from a free fit to the data that C^{ 0 ) is indeed small and € 4 ( 0 ) is close to 
the Adler's value of —0.3 (see Figure p. We use € 4 ( 0 ) = —0.3 in our comparisons 
in the next section. 



CjlO) 


Fig.l. One- and three-standard deviation limits on Cf (0) and 04 ( 0 ) as extracted from 
measurements of "v^p —> The square denotes the model predictions by Adler 

[p^. (Figure adapted after [^.) 


4 Interpretation of Cf (Q^) in the linear a-model 

The axial N-A transition amplitudes can be interpreted in an illustrative way in 
quark models involving chiral fields like the linear cr-model (LSM), which may 
reveal the importance of non-quark degrees of freedom in baryons. Due to dif¬ 
ficulties in consistent incorporation of the pion field, the model predictions for 
these amplitudes are very scarce The present work [ p^JT^ was partly also 
motivated by the experience gained in the successful phenomenological descrip¬ 
tion of the quadrupole electro-excitation of the A within the LSM, in which the 
pion cloud was shown to play a major role 

4.1 Two-radial mode approach 

We have realised that by treating the nucleon and the A in the LSM in a simpler, 
one-radial mode ansatz, the off-diagonal Goldberger-Treiman relation can not be 
satisfied. For the calculation of the amplitudes in the LSM, we have therefore 
used the two-radial mode ansatz for the physical baryon states which allows for 
different pion clouds around the bare baryons. The physical baryons are obtained 
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from the superposition of bare quark cores and coherent states of mesons by the 
Peierls-Yoccoz angular projection. For the nucleon we have the ansatz 

IN) = AAnP^ [®N|Nq) + ONAlAq)] , (7) 


where A/isf is the normalisation factor. Here On and Ona stand for hedgehog co¬ 
herent states describing the pion cloud arormd the bare nucleon and bare A, re¬ 
spectively, and pi is the projection operator on the subspace with isospin and 
angular momentum j. Only one profile for the u field is assumed. For the A we 
assume a slightly different ansatz to ensure the proper asymptotic behaviour. We 
take 


|A) =AtA<^P^OA|Aq) + 


dkq(k)[a^,{k)|N)]H 


( 8 ) 


where A/a is the normalisation factor, |N) is the ground state and [ ] i i denotes 
the pion-nucleon state with isospin | and spin j. We have interpreted the lo¬ 
calised model states as wave-packets with definite linear momentum, as elabo¬ 
rated in [|T^. 


4.2 Calculation of helicity amplitudes 

We use the kinematics and notation of [p^. For the quark contribution to the two 
transverse (A = 1) and longitudinal (A = 0) helicity amplitudes we obtain 
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SaA 


= -(A 
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3 3 1 

„22 
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Here u and v are upper and lower components of Dirac spinors for the nucleon 
and the A, while Cq is a coefficient involving integrals of the function q(k) ap¬ 
pearing in (^. The reduced matrix elements of ctt can be expressed in terms of 
analytic functions with intrinsic numbers of pions as arguments. In all three cases, 
we take sa = f ■ For the scalar amplitudes, we take A = 0 and sa = j, and obtain 


stq) =-(Ai i| dr 


1 1 ; 




drr^ ji (kr) (uaVn - vaUn) (Ab||cTT||N) 


For the non-pole meson contribution to the transverse and longitudinal helic¬ 
ity amplitudes we assume the same cr profiles around the bare states, but different 
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for the physical states. Introducing an "average" cr field cr(r) = -1- crA(r)] 

we obfain 


= (A 


Sa j I 


dre‘'^^ ((cT-fTt)VA7To -tioVact) |Ns^_a^/ 
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dtp AN 2tpAN 

dr r 


do 

^tpAN 


+(3A2-2) 


xCf 


drr^ 12 (hr) 


11 

C 22 


(cj — f. 
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where tpAN = (A | tt | N). To compute the scalar amplitude, we make use of the 
off-diagonal virial relation derived in [13| and define 


We obfain 

S(W = -(A, 
Sti 

^ “T 

By using 


Ti' 


tj’’(r) = 


dre^'^^ ((tj- f„)P„o - P(t7To) INh 
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the quark and non-pole meson contributions to the transverse amplitudes can 
finally be broken into L = 0 and L = 2 pieces. 


( 9 ) 

( 10 ) 

( 11 ) 


Al = y- (A°-A 2 ) , 

ai = Aa| = ^(a»-V), 

^ (A° + 2A2) , 


and inserted into (76), (77), and (78) of 0. The pole part of the meson contribu¬ 
tion is 

[2 
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)(Q")=h 


G7tNA(Q^) 




-6 (pole) W ; 'Ti 2 Mn m2 -I- Q2 V 3 
The strong NA form-factor Gttna can be computed through 


GtinaIQ^) Ma + ATn _ 1 
2 Mn 2 Ma “ Tic ^ 
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where the current J has a component corresponding to the quark source and a 
component originating in the meson self-interaction term (see (58) of 


Jo(r) = jo(r) 


97to(r) 


4.3 Results 

Fig- 1 shows the C^(Q^) amplitude with the quark-meson coupling constant of 
g = 4.3 and ma = 600 MeV compared to the experimentally determined form- 
factors. The figure also shows the Cf^(Q^) calculated from the strong ttNA form- 
factor using the off-diagonal Goldberger-Treiman relation. 



Fig. 2. The amplitude Cf ((^) in the two-radial mode LSM. The experimental imcertainty 
at = 0 is given by Eq. (|^. The error ranges are given by the spread in the axial-mass 
parameter Ma as determined from neutrino scattering experiments (broader range, 0 ) 
and from electro-production of pions (narrower range, Eq. (^). Full curves: calculation 
from helicity amplitudes (^, (^, and (pi|); dashed curves: calculation from Gtina- 


The magnitude of C^(Q^) is overesfimafed in fhe LSM, with C^(0) about 
25 % higher than the experimental average. Still, the Q^-dependence follows fhe 
experimental one very well: the Ma from a dipole fit to our calculated values 
agrees to within a few percenf with the experimental Ma- On the other hand, 
with C^(Q^) determined from the calculated strong ttNA form-factor, the abso¬ 
lute normalisation improves, while the fall-off is sfeeper, with Ma « 0.80 GeV. 
Since the model states are not exact eigenstates of fhe LSM Hamiltonian, the dis¬ 
crepancy between the two calculated values in some sense indicates the quality 






























Axial currents in electro-weak pion production 


75 


of the computational scheme. At where the off-diagonal Goldberger- 

Treiman relation is expected to hold, the discrepancy is 17 %. The disagreement 
between the two approaches can be attributed to an over-estimate of the meson 
strength, a characteristic feature of LSM where only the meson fields bind the 
quarks. 

Essentially the same trend is observed in the "diagonal" case: for the nucleon 
we obtain qa = 1.41. The discrepancy with respect to the experimental value of 
1.27 is commensurate with the disagreement in C^(0). The overestimate of qa 
and Ga(Q^) was shown to persist even if the spurious centre-of-mass motion 
of the nucleon is removed An additional projection onto non-zero linear 
momentum therefore does not appear to be feasible. 

The effect of the meson self-interaction is relatively less pronounced in the 
strong coupling constant (only ~ 20 %) than in Cg (Q^). Both CttnaIO) and GttnnIO) 
are over-estimated in the model by ~ 10%. Still, the ratio G 7 tNA( 0 )/G 7 tNN( 0 ) = 
2.01 is considerably higher than either the familiar SU( 6 ) prediction a/ 72/25 or 
the mass-corrected value of 1.65 and compares reasonably well with the ex¬ 
perimental value of 2.2. This improvement is mostly a consequence of the renor¬ 
malisation of the strong vertices due to pions. 



Fig. 3. The amplitude C 4 (Q^) in the two-radial mode linear u-model, with model param¬ 
eters and experimental uncertainties due to the spread in Ma as in Fig. ^ and in the 
Cloudy-Bag Model (see below for discussion). Experimentally, C 4 (0) = —0.3 db 0.5 (see 
0 and Fig. For orientation, the value for C 4 (0) is used without error-bars. 


The determination of the C 4 () is less reliable because the meson contribu¬ 
tion to the scalar component of this amplitude [ p^ is very sensitive to small vari¬ 
ations of the profiles. However, the experimental value is very uncertain as well. 
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Neglecting the non-pole contribution to the scalar amplitude and Cg (Q^) (with 
the pole contribution canceling out), C 4 (Q^) is fixed to —(M^/2M^) C 5 (Q^). At 
= 0, this is in excellent numerical agreement with (^). In the LSM, the non¬ 
pole contribution to C^(Q^) happens to be non-negligible and tends to increase 
C 4 (Q^) at small Q^, as seen in Fig. ^ An almost identical conclusion regarding 
C 4 (Q^) applies in the case of fhe Cloudy-Bag Model, as shown below. 

The amplitude is governed by the pion pole for small values of and 
hence by the value of Gtcna which is well reproduced in fhe LSM, and under¬ 
estimated by ~ 35 % in the Cloudy-Bag Model. Fig. ^ shows that the non-pole 
contribution becomes relatively more important at larger values of Q^. 



A 0 0.25 0.5 0,75 1 

k=0 Q' 1 (GeV/c)^ 1 


Fig. 4. The non-pole part and the total ainplitude CglQ^) in the two-radial mode linear 
cr-model. Model parameters are as in Fig. 0 


5 Interpretation of Cf (Q^) in the Cloudy-Bag Model 

For the calculation in the Cloudy-Bag Model (CBM) we have assumed the usual 
perturbative form for fhe pion profiles using the experimental masses for fhe nu¬ 
cleon and A. Since fhe pion contribution to the axial current in the CBM has the 
form f 7 t 9 “ 7 t, only fhe quarks confribufe fo fhe C 4 (Q^) and C^(Q^), while C^(Q^) 
is almosf complefely dominated by the pion pole (see contribution by B. Golli 
[[ 1 ^). With respect to the LSM, the sensitivity of fhe axial form-factors to the non¬ 
quark degrees of freedom is therefore almost reversed. 

In the CBM, only the non-pole component of fhe axial currenf contributes 
to the amplitudes, and as a result the C 5 (0) amplitude is less than 2/3 of the 
experimental value. The behaviour of C 5 () (see Fig. ^ is similar as in fhe pure 
MIT Bag Model (to within 10 %), with fitted Ma ~ 1.2 GeV fm/R. The off-diagonal 
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Goldberger-Treiman relation is satisfied in the CBM, but from G^na has a 
steeper fall-off with fitted Ma ~ 0.8 GeV fm/R. 



Fig. 5. The amplitude Cf ( Q ^) in the Cloudy-Bag Model for three values of the bag radius. 
Experimental uncertainties are as in caption to Fig. 


The large discrepancy can be partly attributed to the fact that the CBM pre¬ 
dicts a too low value for Gt^nn/ and consequently Gt^na- We have found that the 
pions increase the Gtina/G^nn ratio by ~ 15% through vertex renormalisation. 
The effect is further enhanced by the mass-correction factor 2 Ma/(Ma + Mn), yet 
suppressed in the kinematical extrapolation of GtinaIQ^) to the SU( 6 ) limit. This 
suppression is weaker at small bag radii R: the ratio drops from 2.05 at R = 0.7 fm 
to 1.60 (below the SU( 6 ) value) at R = 1.3 fm. 

The determination of the C 4 (Q^) is less reliable for very much the same rea¬ 
son as in the LSM. The non-pole contribution to C^(Q^) tends to add to the ex¬ 
cessive strength of C 4 (Q^) at low Q^, as seen in Fig. || Never the less, the exper¬ 
imental data are too coarse to allow for a meaningful comparison to the model. 
For technical details regarding the calculation in the CBM, refer to [p^. 
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